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ABSTRACT 


The  logarithmic  oxidation  equation  y  -  k  ta  (~,+  l)  la  derived 


assuming  control  of  the  rate  by  electron  flow  from  metal  to  oxide.  Electron 


flow  during  oxidation,  in  turn,  is  a  function  of  a  changing  less  positive 
or  a  more  negative  apace  charge  in  the  oxide  extending  up  to  several 


thousand  angstrom  units  from  the  metal  surface,  where  increasing  numbers 


of  trapped  electron!  at  lattice  defect  sites  account  for  the  changing 
apace  charge.  The  apace  charge  is  composed  of  two  parts:  (l)  a  uniform 
charge  density  layer  next  to  the  metal  and  (2)  a  diffuse  charge  density 


layer  beyond  the  uniform  layer.  Oxidation  follows  the  logarithmic  equation 


during  formation  of  both  apace  charge  layers ,  but  not  afterward,  with  a 


higher  oxidation  rate  accompanying  formation  of  the  diffuse  layer.  It  is 


shown  that  under  some  circumstances  the  particular  distribution  of  negative 


charge  in  the  diffuse  layer  may  lead  to  the  cubic  oxidation  equation . 


Experimental  conditions,  especially  inpurities  present  in  the  oxide  and  its 


heat  treatment,  probably  determine  which  electric  charge  distribution  is 


favored.  For  oxidation  beyond  the  space  charge  layer,  either  the  linear 


equation  is  obeyed,  with  control  of  the  rate  still  focussed  at  the  metal* 
oxide  interface,  or,  more  frequently,  the  parabolic  equation  is  obeyed 


with  control  centered  in  diffusion  and  migration  processes  in  the  oxide. 


in  accord  with  Wagner's  theory. 


from  oxidation  data  for  copper,  the;  density  of  trapped  electrons 


I 


in  the  uniform  charge  density  CU2O  layer  at  150s  is  calculated  to  be  1.3  x 
1015  which  decreases  to  6  x  10^  at  250°C.  The  number  of  available  sites 


IP 


for  trapped  electrons  in  the  diffuse  layer  is  calculated  correspondingly  to 
be  1.0  x  10^  and  2.6  x  10^  respectively.  The  thickness  of  the  uniform 
i&arge  density  layer  varies  from  350  A  at  150s  to  1630  A  at  250*C. 


',  J  The  Rideal- Jones  empirical  relation  AB  *  $  -  3.6,  where  AE  is  the 
activation  energy  for  oxidation  and  0  the  metal  work  function,  both  in 
electron  volts,  and  which  has  been  shown  to  hold  for  carbon,  tungsten  and 
platinum  is  shown  herewith  to  apply  to  nickel,  tantalum  and  titanium,  and 
more  approximately  to  copper  and  iron.  This  relation  is  derived  theoretically 
from  the  same  fundamental  assumptions  used  in  deriving  the  logarithmic 

1 

oxidation  equation. 
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Whenever  the  oxidation  rate  of  a  metal  is  controlled  by  diffusion  of 
Ions  (or  lattice  vacancies)  and  by  migration  of  electrons  (or  positive  holes) 
through  the  oxide  film,  the  familiar  parabolic  aquation  applies s 

y*  -  kxt  +  C  (1) 

where  y  la  the  thickness  of  the  film,  t  la  time,  and  ki  and  C  are  constants, 
(C  may  sometimes  be  aero.)  Wagner established  the  fundamental  reaction 


C.  Wagner,  Z.  physik. 


21B.  25  (1933);  52,  447  (1936);  4GB,  455  (1938). 


kinetics  for  this  equation  and  was  able  to  calculate  the  rate  constant  k, 
and  its  dependence  on  caqFgsn  partial  pressure,  from  independent  physical 
chemical  properties  of  oxide  and  metal. 

However,  oxidation  In  the  lower  temperature  range  or  when  thin  films 
are  formed  Is  found  to  follow  the  so-called  logarithmic  equation: 

7  *  kQ  A  (£  +  1)  (2] 

where  k  and  are  constants.  For  t))?’,  this  equation  is  approximated 


and  y  becomes  linear  plotted  with  logarithm  of  t. 

then  the  logarithmic  aquation  was  first  announced  in  19&  by 


to) 

and  KfJster'  '  on  empirical  ground* ,  It*  validity  was  quest! oced.  However, 

Cp  A . . .  -  -  ■  ■  . . — - . ; — -  1  ■  - - - 

'  y0.  Tamoann  and  W.  Kflster,  Z.  Anorg.  Allg.  Chem. ,  123.  196(1922). 

careful  measurements  since  then  have  established  beyond  doubt  that  for 
oxidation  of  many  metals  under  specific  conditions  of  temperature  and  time, 
this  equation  and  no  other  describes  the  facts.  The  logarithmic  equation, 
for  example,  has  been  shown  to  apply  to  oxidation  of  copper^ 2»3*4>5>6)^ 

^B.  Luatman  and  R.  Mchl,  Trans.  AIMS.  143.  246  (l94l). 

(h)A.  H.  White  and  L.  Qermer,  Trans.  Electrochem.  Soc.,  81,  305  (1942). 

^H.  H.  Uhlig  and  A.  Brenner,  Acta  Metallurgica,  2 >  10B  (1955). 

^6^R.  lyiecote,  J.  Inst.  Metals,  38,  327  (1950). 

zinc^2|^\  cadmium^2),  lead^,  tin^2^,  manganese^2),  iran^2'®'?*10'11^, 

'^W.  Vernon,  *.  Akeroyd  and  B.  Stroud,  J.  Inst.  Metals,  6^,  301  (1939). 

^B.  Luatman,  Trans.  Electrochem.  Soc.,  81,  359  (1942). 

^A.  Winterbottom,  J.  Iron  and  Steel  Inst.,  165.  9  (1950). 

^10^W.  Vernon,  E.  Calnan,  C.  Clews  and  T.  Nurse,  Proc.  Royal  Soc.,  2l6A,  375 
(1953). 

Davies,  U.  Evans  and  J.  Agar,  Proc.  Royal  Soc.,  225A.  443  (1954). 
nickel^2'12'15),  aluminum^5'1^,  titanium^,  and  tantalum^16). 

(12V  H.  Uhlig  and  J.  J.  Pickett,  unpublished  data.  . 

(15\f.  Scheuble,  Z.  Fhysik.,  1^,  125  (1953). 

Steinheil,  Ann.  Fhysik.,  lg,  455  (1934). 

(15^J.  Waber,  0.  Sturdy  and  E.  Wise,  22&  (1953) . 

Waber,  0.  Sturdy,  E.  Wise  and  C.  Tipton,  J.  Electrochem.  Soc.,  gg,  121  (1952). 


•  2  • 


Several  attempt!  have  been  made  to  describe  a  mechanism  of  diffusion  through 

oxides  resulting  in  the  logarithmic  equation,  such  as  by  assumed  specialised 

(17) 

flaw-paths  or  sones  of  loose  structure  in  oxides''  1 7,  or  by  diffusion  blocks 

WE  r.  Evans,  Trans.  Electrochem.  Soc.,  8%,  335  (1953). 

and  leakage  paths  in  the  oxide (**#*8)^  or  by  an  assumed  Ion  concentration 


™U.  B.  Evans,  Nature.  157.  732  (1946). 


gradient  or  polarization^  Mott ^ 20 ^  contributed  an  analysis  of 


Miley,  Trans.  Electrochem.  Soc.,  81.  391  (1952).  . 

0.  Hamilton  and-H.  Miley,  Ibid,  §1/  513  (1952). 

^U.  Mott,  Trans.  Faraday  Soc.,  21>  1175  (1939). 

♦Vernon:  Chews' ^ 

as  found  to  apply  to  thin  film  .oxidation  of  iron,  expresses  oxygen  diffu¬ 
sion  through  the  oxide  to  the  metal  surface,  whereas  the  parabolic  equation, 
‘when  it  applies,  corresponds  to  diffusion  of  iron  ions  through  the  oxide 
from  metal  to  gas.  They  base  their  opinion  on  observed  greater  weight  gain 
of  iron  specimens  when  oxidized  in  the  logarithmic  rate  region  then  corres¬ 
ponds  to  oxygen  in  the  Isolated  oxide  filnv  in  contrast  to  correspondence 
of  weight  gain  and  oxygen  in  the  oxide  in  the  parabolic  region.  Their 
observations,  however,  can  be  explained  by  surface  oxidation  of  carbon 
contained  in  their  iron  (0.09^  C),  occurring  parallel  with  the  oxidation 
of  iron,  the  carbon  oxidation  products  being  retained  by  the  iron  oxides 
perhaps  aa  carbonates.  For  higher  temperatures  and  thicker  exile  films 
where  the  parabolic  equation  holds,  carbonates  are  not  stable  and,  hence, 
are  not  expected  to  be  retained  by  the  film.  Oxidation  of  carbon  also 
explains  the  delayed  appearance  of  interference  colors  which  they  report. 

In  this  connection,  the  effect  of  carbon  probably  enters  as  a  factor  in 
the  thinning  (by  reduction)  of  oxide  films  on  iron  during  vacuum  treat¬ 
ment  reported  by  Davies,  Evans  and  Agar(ll).  A  pertinent  reference  in 

this  regard  is  the  work  of  V.  Hola^21^. 


Holm.  J.  Res.  Eat.  Bur.  Side..  26,  569  (1958). 
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“*  10°tC  ftmU"  *  «*■  *“*•««.  fU»  of  etoinhio- 

“trt0  --«.«*  w  .  *»***.  .^tl0B.  „ 

y**  m  *“ 1  *"**  -  «—  — to  .*P"..  rat.  Of 

aiT  rr *•*■“• ,**  *•  “*”•*  c°otrou*4 *•  -«*«- 

PM  *0,1  accounted,  houever,  for  oxide  file.  ally  in  the  coder  of 

'  *•  r  ^  »-**•  T—  *•  P-tlo.  hold.  f.r  fu. 

“  *“*  “  10’000  A-  -  —  *•  Mta  VMM*  a  t-perature  lade- 
pendent  ko,  contrary  to  shat  in  observed.  Subeegiently,  Nett  and  C.brer.f®*8’) 

^N'  Mott»  Tr“8‘  Farad»y  floe.,  429  (19*7)1  - - - - 

(23) 

_  I*t*  ******  *p°rte  on  Prog,  in  Ayeics,  I2j  165  (19119). 


hypothesized  an  alternative  rat.  octroi  *  diffusion  of  act*  iM  ^ 
than  escape  of  electron.,  assualng  that  negatively 
adsorbed  on  tta  oxide  surface  create  an  electric  field  wthln  th.  oxide 
^loh  induce,  algratlon  of  PO.itlv.  l»s  to  tb.  oxide  eurf„e.  ^  ^ 
that  under  condition.  *ere  IgTstioa  velocity  of  pontiv.  1=0.  1. 
tloul  to  field  etrength,  a  p^sbollc  oadatlca  entice  r.ente,  but  far 
thinner  files  far  idiicb  the  atgratlea  velocity  le  proportional  to  an 
exponential  toner  of  the  field  strength  end  rate  of  ..cape  of  eetal  lone 
Into  the  oxide  beooaea  ocotmllla*  an  e*»tloo  oftt.twe£-0-hlatl. 
henwd  (the  Inverse  logaritbalc  equation),  ai.  .1.0  led  to  the  ,ltu.tloa 

*”  *"»”’**'“>  *  «-«*>.  arid,  file  thlctn...  could 

ho  expected.  At  the  legndanc  e*.tl«  (*>  ob«r^  ,*«**«* 

not  reault  from  thaae  partictilar  M«u^>tloDS. 

Caapbell  end  IhcaM^^)  erpreeeed  the  opinion  that  laltui  deviation. 

C.apb.11  and  0.  Ihoeae,  Ben,.  Ilectroohea.  doc.,  Jlj.  &3  (19I17). 
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from  the  parabolic  equation  of  oxidation  for  copper  are  caused  by  a  space 

4  0 

charge  layer  in  the  oxide  about  10  A  thick  as  determined  by  electrical 

capacitance.  Tylecote^  confirmed  that  the  logarithmic  equation  holds  for 

5  ° 

films  less  than  about  10  A  thick,  whereas  the  parabolic  equation  applies 

t>U  sm**^ 

to  thicker  films.  Later  discussion  in  taaa  paper  bears  out  that  the 
logarithmic  equation  is  related  to  conditions  applying  famSmg  formation  of 
the  space  charge  layer,  and  that  the  mechanism  of  oxidation  changes  for  thicker 
oxide  films. 

Oxidation  Control  at  the  Metal-Oxide  Interface 

It  is  perhaps  obvious  that  an  acceptable  mechanism  of  thin  film  oxida¬ 
tion  of  metals  should  first  establish  the  actual  rate  controlling  process. 

A  detailed  examination  of  reaction-rate  data  provides  evidence  that  the 
rate  is  controlled  largely  by  processes  at  the  metal -oxide  Interface  rather 
than  by  reactions  at  the  oxide -oxygen  interface,  or  by  diffusion  through 
the  oxide.  This  evidence  comes  from  the  effect  on  the  oxidation  rate  oft 
(l)  crystal  orientation,  one  crystal  face  oxidiring  at  a 

rate  appreciably  different  from  any  other  crystal  face^®*2^. 

^^A.  Qvathmey,  "Corrosion  Handbook",  p.  33-55#  edited  by  H.  E.  fihlig, 

John  Wiley  and  Sons,  Inc.,  Rev  York,  1946. 

Luatman^,and  Me  hi  and  McCandlees^2^  pointed  out  that 


(26Jr.  Mehl  and  E.  MeCendless.  Trans.  A.I.M.E.,  125.  531  (1937). 

these  differences  cannot  be  ascribed  to  differing  diffusion 
’  rates  through  oxides  specifically  oriented  on  metal  crystal 
faces,  because  the  oxides -of  iron  and  copper  at  low  tempera¬ 
tures  ere  cubic,  and  diffusion  through  cubic  lattices  is 


Isotropic.  Differing  degree  of  sintering  or  differing  re* 
crystallization  of  oxide  as  a  function  of  crystal  face  is 


also  not  an  acceptable  explanation,  because  Mehl,  McCandless 
(27) 

and  Shines'  ' '  found  that  oxides  grown  on  single  crystals  of 
^R.  Mehl,  1.  McCandless  and  F.  Shines,  Nature,  134,  1009  (1934). 


copper  remain  single  crystals  even  after  they  have  developed 
into  thick  scales,  precluding  the  possibility  of  any  sinter* 
ing  or  recrystallization. 

(2)  Lattice  Transformation 

The  activation  energy  for  oxidation  of  iron  In  CO®  + . 
HgO  undergoes  a  discontinuity  at  the  transformation  tempera¬ 
ture  of  iron  from  a  (bcc)  to  7  (fee)  lattice A  similar 

^K.  Fischbeck  and  F.  Salzer.  Metallwirtschaft ,  l4,  753  (1935). 

(29) 

discontinuity  vas  found  for  a  IO.65&  chromiuin-iron  alloy ' 


*H.  H.  Uhlig  and  A.  8.  Brasunas,  J.  Electrochem.  Soc.,  gj,  448  (1950) 


heated  in  oxygen  under  conditions  where  thin  films  are  formed 


(3)  Curie  Temperature 

The  activation  energies  for  oxidation  of  several 
chromium-iron  alloys  in  oxygen  are  higher  above  the  Curie 
temperatures  than  belov^2”^.  This  relation  is  also  found 
to  hold  for  nickel^12).  Taranann  and  Siebel^^  furthermore 


Wo.  Thomson  and  0.  Siebel,  Z.  anorg.  allg.  Chen.,  148,  297  (1925)* 

.  -  ,  .  —  - - -  —  -  ,  1 

found  discontinuities  in  the  oxidation  rates  of  iron-nickel 
alloys  at  tbs  Curie  temperatures.  All  these  observations 


show  that  magnetic  transitions  in  astals  at  the  Curie  tenpers- 
tttra  can  affect  the  oxidation  rate. 

The  overall  evidence  for  control  of  rate  at  the  netal- 
cadde  interface  applies  whenever  thin  films  are  toned.  For 
thicker  films,  of  course,  where  diffusion  through  the  oxide 
is  contro] .  'ng,  factors  of  grain  orientation,  lattice  trsna- 
formation  and  magnetic  chugs  are  mo  laager  lsvoptant. 

Control  by  Electron  Balssion 

110^(22,23)  favored  slow  escape  of  metal  ions  from  the  metal  into  the 
oxide  lattice  on  the  controlling  step  in  several  reactions  producing 
relatively  thin  films,  and  electron  ealsslon  as  controlling  in  insulator 
type  films  belov  kO  A  thick.  Convincing  evidence  is  available  to  show 
that  actually  electron  flow  from  metal  to  ood.de  is  the  alow  step  for  several 
metals,  and  therefore  dominates  the  reaction  process  in  formation  of  films 
as  thick  as  several  thousand  Angstroms.  This  evidence  in  part  cosies  from 
the  empirical  relation  first  proposed  by  Rldeal  and  Hsnsbrot4$~Jbnes^'1'). 

^^B.  Rldeal  and  0.  Waasbrougb- Jones .  Proc.  Royal  Boc. ,  123A.  202  (1929). 

They  showed  that  for  oxidation  of  platinum,  carbon  and  tungsten  the  follow¬ 
ing  relation  holds: 

dl  «  0O  -  K  (k) 

whore  AB  is  the  observed  activation  energy  for  oxidation,  and  0  is  the  wqrit 

w  • 

function  of  tarn  metal  or  carbon.  The  term  K  is  a  constant  equal  to  3.6 
electron  volts  (83  kcal),  and  was  expressed  by  Rldeal  and  Wansbrough -Jones 
as  the  sum  of  4  where  is  the  energy  of  adsorption  and  is  the 
|l|ctroo  affinity  of  the  adsorbed  oxygen  molecule.  In  other  words,  tor 
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these  elements,  each  of  whose  oxides  .at  elevated  temperatures  is  volatile, 
the  activation  energy  for  oxidation  is  related  directly  to  ties  energy 
necessary  to  transfer  an  electron  from  metal  to  oxygen. 

According  to  Iq.  (4),  metals  with  work  function  less  than  3*6  electron, 
volts  (e.g.,  Ha  and  Ca)  possess  an  activation  energy  for  oxidation  equal  to 
zero  (it  cannot  he  negative),  and  the  metal  oxidizes  initially  as  rapidly 
as  oxygen  comes  into  contact.  For  metals  of  higher  work  function,  an  the 
other  hand,  (e.g.,  Hi,  Pt,  W),  the  activation  energy  is  positive,  the  reac¬ 
tion  is  delayed,  and  more  oxygen  reaches  the  metal  surface  than  can  react. 

This  means  that  oxygen  is  able  to  adsorb  (chemisorb)  on  the  meted,  for  which 
it  has  a  certain  affinity,  remaining  there  a  measurable  time  before  metal 
atoms  leave  their  lattice  to  Initiate  an  oxide  lattice. 

It  is  of  Intel*  v  to  note  that  the  measured  value  of  K  in  Xq.  (4-) 

(3.6  e.v. )  approximates  the  electron  affinity  of  the  oxygen  atom  (3>fc  e.v.)^2\ 

"Oxidation  Potentials",  W.  M.  latlmer,  second  edition,  p.  lB,  prentice  - 
Hall,  Inc.,  Sew  York,  1952. 

The  small  difference  of  about  0.6  electron  volt  (4600  cal.)  ia  possibly 
accounted  for  by  the  sum  of  energy  changes  associated  with  adsorption  and 
dissociation  of  molecular  oxygen  on  the  metal  surface.  In  other  words,  the 
alov  process  of  oxidation  of  carbon,  Pt  and  W  can  be  Interpreted  as  the 
transformation  of  physically  adsorbed  molecular  oxygen,  in  which  only  meek 

4 

% 

Van  der  Waals  binding  forces  apply,  to  chemisorbed  atomic  oxygen  with  the 
meted  acting  aa  electron  donor  end  oxygen  as  electron  acceptor.  Accordingly, 
equation  (4)  indicates  that  both  rates  of  physical  adsorption  of  oxygen 
on  the  metal  surface  and  volatilization  of  the  metal-oxygen  ion  complex,  at 

A  *  k  . ij- ,  *  .»  ‘  •  ,  , 

the  high  temperatures  presently  considered,  are  rapid  compared  to  uheml.  sorption 
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of  oxygen .  This  conclusion  Is  in  agreement  with  previous  expressions  to  the 


effect  that  chemisorption  is  often  (hut  not  always)  a  slow  process  requiring 
an  activation  energy^*  It  also  agrees  with  previous  proposals  that 

f"t«V . .  ■  . . . : . . . .  — ■  . . — 

S.  Taylor,  J.A.C.8.,  ^  578  (1931). 

^^"Chemisorption",  B.  Trapcell,  Academic  Press,  Hew  York,  1955* 

the  process  of  physical  adsorption  precedes  eheai sorption ( 3 5 » 36 ) ^ 

Becker  and  C.  Hartman,  J.  Fhys.  Chem.,  57.  153  (1953). 

^56^G.  Hirlich,  J.  ays.  Chem.,  473  (1955 )• 

I  .  - 

. .  mine  is  ■  — fcm  m  m  ■  ■...!■  ■■■■■—  .  . . .  mmmUmmmmrnm e— mm . .  ■ 

I  * 

It  is  significant  that  the  Ridsal-Jtaaes  relation  eaa  he  applied  also 
to  some  metals  having  non-volatile  oxides  under  conditions  where  the  oxide 
films  are  thin  and  the  logarithmic  equation  is  obeyed.  Hie  data  are 
summarised  In  Table  I,  and  provide  evidence  that  for  such  metals  oxidation 
continues  to  be  controlled  by  electron  flow  from  the  metal.  Presumably, 
the  energy  gained  by  the  electron  entering  the  oxide  is  almost  exactly 
compensated  for  by  the  electron  leaving  the  oxide  again  at  the  oxide -air 
Interface  to  enter  the  oxygen  atom,  now  adsorbed  not  on  the  metal,  but  on 
the  oxide.  The  negatively  charged  oxygen  ions  eventually  enter  lattice 
positions  in  the  oxide  in  accord  with  the  Wagner  mechanism  of  oxidation 
and  tarnish.  Therefore,  AE  is  again  given  essentially  by  the  difference 
of  the  metal  work  function  and  the  electron  affinity  of  oxygen.  From  the 
above  considerations,  Iq.  (4)  can  be  expected  to  hold  whenever  electron 
escape  from  the  metal  controls,  and  when  the  work  function  of  the  oxide 
has  approximately  the  earns  value  at  the  metal -oxide  and  oxide -environment 
interfaces. 

The  observed  relation  between  work  function  and  activation  energy 
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TABLE  I 


hslakqr  between  activation  ENERGY  FOR  OXIDATION 
_ AND  WORK  FUNCTION:  &  -  0  - 


£Z 

0 

ro 

K 

Pt 

w 

c 

2.74  e.v. 
O.87 

0.52 

6.35  e.v.* 
4.48* 

4.31* 

3.61  e.v. 
,.t  3-61 

3.79 

m 

Ta 

H 

0.90» 

0.55a 

0.56o 

4.50** 

4.12*** 

3*95** 

3.60  -x.-x- 

3.57 

3.39 

WAese  values  employed  by  Rl&al  and  Wansbreugh- Jones  are  close  to  mean  ‘ 
valixes  given  by  H.  Michael  eon  baaed  on  a  review  of  the  literature,  with 
the  exception  of  platinum  where  the  value  in  Table  I,  alloying  photo¬ 
electric  data  of  DuBzldge  (1928)  is  1.06  volt  higher! 

B  H*  H*  Uhll8  and  J.  J.  Pickett,  unpublished  data, 
a  Beference  16. 
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c lari fie*  vtay  the  rate  of  oxidation  should  vary  vlth  crystal  faoe,  since 
the  work  function  also  varies  with  crystal  Furthermore, 

^R.  Snoluchovricl,  ftys.  Rev.,  60,  66l  (19*H). 

Farnsworth  and  R.  Winch,  Hiys.  Rev.,  %8,  812  (1940). 

^^Also  various  papers  on  field  emission  microscope,  e.g.,  I,  W.  MtiUer, 
Ergeb.  der  Exakten  Naturwias . ,  2JJ  290-360  (1953/ J  J*  A.  Becker,  Bell 
Syst.  Tech.  J.,  2£,  907  (1951). 


the  lank  function  of  7  Iron  differs  from  that  of  a  iron.  In  accord  vith 

change  of  activation  energy  for  oxidation  above  and  below  the  transition 

(ko\ 

temperature.  Boudremont  and  Rfidiger'  '  reported  a  higher  work  function 

Boudremont  and  0.  Hfldlser.  Ifaturwles.,  39.  399  (1952). 
for  7  than  for  a  Iron,  but  Wahlin^)  reported  the  reverse,  namely,  a  higher 

Wshlln,  aye  Rev.,  &,  509  (19^2). 
value  for  a  than  for  7  iron*. 

Along  the  seam  lines,  values  of  vork  function  for  nickel  above  the 
Curie  tesgperature  (350°C)  are  higher  than  values  below  the  Curie  temperature 

Cardwell.  Ron.  Rev..  76,  125  (1949). 

corresponding  to  an  observed  activation  energy  for  oxidation  of  nickel  In 
oxygen  above  the  Curie  teaperature  that  la  1150  calorie s/mole,  or  0.05  e.v. 
greater  than  the  value  below  the  Curie  temperature  These  correlations 


point  toward  the  iaportent  part  played  by  electron  flov  from  the  netal  la 
the  Initial  oxidation  process,  and  focus  the  necessity  of  attention  on 
various  f actors  affecting  electron  emission  during  oxide  growth* 

Space  Charge  Effects 

If  at  equilibrium,  an  oxide  has  greater  electron  affinity  or  higher 
work  function  than  the  aetal  In  contact  from  which  it  forms,  the  oxide  will 
tend  to  acquire  excess  negative  charge,  and  the  metal.  In  turn,  will  acquire 
an  equal  positive  charge,  resulting  In  an  electrical  double  layer.  Similarly, 
If  the  electron  affinities  are  reversed,  the  oxide  will  lose  electrons  to 
the  metal.  Hie  amount  of  charge  eventually  transferred  across  the  boundary 
In  either  event  will  set  up  a  field  compensating  exactly  for  the  differences 
In  electron  affinity  of  the  two  phases.  Or  stated  another  way,  when  contact 
Is  made  between  a  metal  and  a  seal -conductor,  all  the  electronic  energy 
levels  of  the  eeml -conductor  are  altered  relative  to  those  of  the  metal 
by  the  amount  of  the  contact  potential  difference. 

For  two  dissimilar  metals  in  contact,  the  electric  charge  constituting 
the  double  layer  la  confined  to  a  very  small  region  of  atomic  dimensions  at 
the  Interface.  In  semi -conducting  oxides,  on  the  other  hand,  with  filled 
energy  levels  and  relatively  few  electrons  excited  into  the  conduction  band, 
exeees  charge  associates  Itself  mainly  at  lattice  imperfections,  such  as 
occur  at  Impurity  centers  or  at  lattice  vacancies.  Negative  charge  of  this 
kind  exists  as  so-salbsd  bound  or  ’trapped"  electrons^).  Similarly, 

Mott  and  B.  Chutney,  "Electronic  Processes  in  Ionic  Crystals",  pp.  80- 
88,  Oxford  Press  (19*0). 

bound  or  trapped  slsotrcn  <  deficient  sites  or  positive  holes  asy.Ieffi#t«  X» 
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view  of  the  restricted  ouaber  of  sites  for  trapped  electrons  or  positive 
holes,  the  spsee  charge  la  oxides  extends  over  a  greater  voltaee  than  in  the 
ease  of  metals.  Volta  or  contact  potential  measurements  for  oxidised  copper 
end  sine,  as  described  later,  indicate  that  the  space  charge  nay  extend 
into  the  oxide  film  for  several  thousand  Angstroms*  Vjben  the  oxidation 
rate,  therefore,  is  controlled  by  electron  flow  from  the  metal,  an 
increasing  negative  space  charge  (or  decreasing  positive  space  charge) 
obviously  enters  as  an  important  factor  in  establishing  the  rate  of 
electron  flow,  tending  in  general  to  slow  down  escape  of  electron*  as  the 
oxide  grows,  and,  therefore,  simultaneously  slowing  deem  escape  of  metal 
ions  and  the  oxidation  process  itself. 

Volta  Potential  Measurements  and  the  Initial  Oxidation  Reaction 

The  oxidation  process,  of  course,  is  not  an  equilibrium  situation, 
but  rather  a  dynamic  building  up  of  oxide  through  interaction  of  the  metal 
and  oxide  with  molecular  oxygen.  The  various  steps  in  the  initial  process 
appear  to  bet 

(1)  rapid  physical  adsorption  of  molecular  oxygen  on  the  metal 

(2)  dissociation  of  moleaular  oxygen  into  shamteoaHad  atomic 

I 

oacygtn 

(£}  sublimation  of  the  metal-oxygen  oosgtlex  to  initiate  an  oxide 
lattice 

(h)  vaporisation  of  metal  lomt  gad  electrons  into  the  oxide  to 
form  additional  oxide 

(5)  formation  of  lattice  vacancies  by  oxygen,  or  combination  of 
interstitial  metal  ions  with  oxygen  at  the  oxide -air  inter¬ 
face.  (Vfe  not  treat  the  relatively  few  known  oases 

where  oxygen  diffuses  to  the  metal -oxide  interface.) 
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If  step  4  controls,  In  accord  with  present  considerations,  the  rate  Of 

•W 

escape  of  ions  Is  given  by  the  expression  K'  exp~  ,  where  W  is  the  heat 
of  sublimation  of  the  ions  from  the  metal  into  the  oxide  at  the  metal-oxide 
Interface,  k  is  Boltzmann's  constant,  T  is  the  absolute  temperature,  and 
K'  is  a  constant.  Electrons  can  be  assumed  to  evaporate  simultaneously 
at  a  rate  equal  to  K"  exp  ,  vhere  0  is  the  vork  function  of  the  metal 
as  modified  by  the  oxide  in  contact.  When  electron  emission  is  the 
slower  process  and  therefore  the  controlling  factor  in  the  formation  of 
oxide,  as  has  been  discussed,  and  since  the  exponential  terms  and  not  K* 
and  K"  in  the  above  rate  expressions  are  usually  dominant,  It  follows  that 
0  in  this  instance  must  be  greater  than  W. 

This  state  of  affairs  leads  to  an  initial  positive  space  charge  In 
the  oxide  when  a  metal  Is  first  oxidized  because  of  the  rapid  escape  of 
positive  metal  ions.  Their  rate  of  escape  almost  immediately  is  established 
by  the  slower  electron  flow  from  metal  to  oxide.  As  the  oxide  grows,  a 
decrease  of  positive  space  charge  occurs  because  a  portion  of  the  electrons 
escaping  Into  the  oxide  become  trapped  at  oxide  lattice  Imperfections . 
Eventually,  the  oxide  may  acquire  an  overall  negative  Charge  depending  on 
the  final  equilibrium  state  of  electric  charge. 

As  oxidation  continues,  rate  of  escape  of  ions  and  of  electrons, 
approach  each  other,  and  some  factor  other  than  electron  flow  assumes 
control  of  the  rate,  particularly  after  the  oxide  grows  to  a  thickness 
corresponding  to  equilibrium  electric  charge  distribution  at  the  metal -oxide 
interface  (or  when  the  Fermi  levels  in  metal  and  oxide  become  identical). 

Thicker  oxides  forming  beyond  the  boundaries  of  the  space  chargs  layer 

.  *  . 

are  electrically  neutral.  If  escape  of  electrons  and  leu  should  continue 

•  1^  • 


- - ' _ - _ _ 

to  control  the  rate  at  this  stage,  as  is  possible,  the  oxidation  rate  would, 
bacons  constant,  leading  to  the  so-called  linear  equation,  (y  »  k*t  ♦  03}* 

But  the  constant  oxidation  rate  In  this  Instance  would  be  preceded  by  an 

(k6) 

initially  more  rapid  and  gradually  decreasing  rate.  Oulbransen  and  tysong'  ' 
^1.  Oulbreneen  and  W.  tfrsong,  J.  ihya.  Chen.,  51,  100?  (19*7). 

i,  ■«  n  . . mmmmmtmmmrnm —  ■  ».—■■■■,  — aee—  n  mm e— me—— — . . 

found  aluminum  to  behave  somewhat  along  these  lines  which  any  confix*  that 
this  type  of  oxidation  control  is  sonatinas  found.  Afanrinua  whan  oxidised 
in  low  pressure  oxygen  below  *50*C  followed  the  parabolic  equation,  and 
above  *50*  C  followed  the  linear  equation.  Whether  one  equation  or  another 
is  obeyed  probably  depends  on  the  defect  structure  of  the  oxide  under  the 
particular  conditions  of  any  given  experiment.  That  Is,  factors  be  cons 
Important  which  determine  whether  site#  continue  to  be  available,  as  the 
oxide  grows,  for  trapping  of  electrons  or  positive  holes,  and  which  establish 
whether  control  of  the  rate  remains  at  the  natal -oxide  interface  or  la 
transferred  to  diffusion  and  conduction  processes  in  the  oxide.  Also, 
patter  is  of  crack  formation  in  the  oxide  aay  play  a  role  in  accord  with  the 
conventional  explanation  of  mechanism  applied  to  the  linear  equation. 
Nevertheless,  It  Is  apparent  that  the  linear  equation,  as  described  above, 
may  also  hold  for  metals  covered  by  a  continuous  protective  oxide  film 
free  of  cracks.  Such  protective  films  have,  in  fact,  been  observed  whan 
nagnealun,  which  follows  the  linear  aquation,  oxidises  initially^***)* 

Ieontis  sad  1.  mines,  Trans.APg,  166.  263  (19*6). 

Oulbransen,  Trans.  KUetrochen.  80c.,  §£>  509  (l9*5). 

Also,  tfcs  black  cmstaosoos  NpO  found  by  Ieontis  and  Mttass^  ant  bo  the 
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metal  surface,  over  which  white  HgO  forms,  may  be  visible  evidence  of  the 
space  charge  layer,  as  veil  as  evidence,  as  they  suggest,  of  a  nou- 
•toichlometric  composition. 

Evidence  for  the  trend  of  space  charge  with  oxide  growth  is  provided 
by  contact  or  Volta  potential  measurements  of  metals  undergoing  oxidation 
in  air  at  elevated  temperatures,  typical  potentials  of  copper  oxidised  in 
air  at  225*C  and  of  sine  at  200*C  are  given  in  Figs.  1  and  2.  these 
measurements  were  made  employing  the  vibrating  condenser  method^*^.  A 


'V.  Zisman,  Rev.  Bel.  Instr.,  367  (1932). 

^H.  H.  Uhllg,  J.  Applied  ttys.,  22,  1399  (1951). 


silver  electrode  as  reference  was  heated  electrically  to  the  same,  temperature 
as  the  copper  or  zinc  electrodes,  all  of  the  same  size  (l  x  1"),  the  tempera* 
ture  of  the  electrodes  being  measured  and  controlled  automatically  by 
means  Of  embedded  thermocouples,  the  initial  oxidation  process  for  copper 
results  in  a  Volta  potential  0.5-0. 7  volt  more  active  than  silver,  and  for 
zinc  about  1.5  volt  more  active  than  silver,,  corresponding  to  buildup  of 
a  positive  space  charge  in  the  metal  oxides.  As  the  metals  oxidize 
further,  these  potentials  become  first  rapidly  more  noble  (more  like  Ag), 
indicating  that  the  positive  space  charge  is  being  neutralized,  followed 
by  a  less  rapid  ant  prolonged  trend  in  the  noble  direction. 

Similar  conclusions  may  be  derived  from  Volta  potential  measurement! 
reported  by  Hirschberg  and  Lange They  tested  mine  surfaces  between 


^  Mrschberg  andLX.  Dangs,.  Xhturwlssenschaften,  Jftl  (1952). 
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Fig-  X  Volta  Potential*  of  Caspar  OKidlslaf  in  Air 
Silver,  lath  at  225 *C. 


CO  N*  iD 


6  6  6  6 

**|oa 


0.4 


! 


* 


20  and  1*07*0  for  20  minutes  at  each  temperature,  and  measured  Volta  potentials 
of  such  surfaces  (presumably  at  room  temperature  rather  than  at  elevated 
temperature  as  above)  with  respect  to  zinc.  At  lover  temperatures  of 
oxidation  and  hence  for  the  thinner  films,  the  sign  of  Volta  potential 
corresponded  to  a  positive  space  charge  in  the  oxide,  Hat  for  higher  tempera¬ 
ture!  and  thicker  oxide  films  the  sign  reversed. 

Volta  potential  measurements  of  metals  undergoing  oxidation  are  not 
easily  reproducible,  so  that  quantitative  features  of  Figs.  1  and  2 
should  not  asnoms  great  weight.  Qualitatively,  however,  the  trend  of 
effects  can  be  reproduced  from  one  run  to  the  next.  These  measurements 
are  being  continued  at  various  temperatures  and  for  various  metals  and  will 
be  reported  elsevhere .  In  all  instances  observed  so  far,  the  sign  of 
Volta  potential  after  an  initial  period  of  oxidation,  corresponds  to  decreas¬ 
ing  positive  space  charge  or  increasing  negative  charge.  Interpretation 
of  Volta  potential  measurements  as  described  above  assumes  that  electric 
charge  at  the  air-oxide  interface  remains  essentially  constant,  idiich  is 
a  reasonable  assumption  for  constant  partial  pressure  of  oxygen;  and  that 
major  change  of  potential  is  accounted  for  by  changing  electric  charge  in 
the  thin  oxide  films. 

The  potential  at  the  metal-oxide  interface  set  up  by  changing  space 
charge  in  the  oxide  layer  alters,  of  course,  the  work  necessary  to  carry 
a  positive  ion  or  electron  across  the  interfaoe,  and,  hence,  influences 
their  rate  of  escape  into  the  oxide.  It  is  this  sgpeeesh  to  the  oxidation 
process  that  is  explored  next. 

Calculation  of  Oxidation  Rate 

Mb  consider  first  that  excess  charge  in  the  growing  oxide  film 
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of  total  thickness  y  is  uniformly  distributed  throughout  the  oxide  volume. 
This  corresponds  to  an  oxide  vlth  uniform  defect  or  impurity  density 
at  which  electrons  or  positive  holes  are  bound.  When  a  film  of  this  kind 
readies  a  supposed  thickness  {,  wb  shall  further  assume  that  the  total 
negative  charge  in  the  film  creates  a  field  at  the  metal  surface  compensat¬ 
ing  exactly  for  the  positive  field.  Hie  positive  field  results  from  eycess 
metal  ions  of  positive  charge  escaping  into  the  oxide,  or,  alternatively, 
results  from  the  electron  "sink"  at  lattice  vacancies  formed  at  the  outer 
oxide  surface.  Therefore,  when  the  hypothetical  film  thickness  becomes 
the  electric  field  at  the  metal  surface  reaches  the  equilibrium  value  of 
sero  (Feral  levels  in  oxide  and  metal  are  equal),  corresponding  to  an  order 
of  film  thickness  beyond  which  the  space  charge  can  exert  either  a  constant 
or  a  negligible  effect  on  electron  escape  and  rate  of  oxide  growth..  If 
n  is  the  density  of  trapped  electrons  in  the  oxide,  and  -e  is  the  charge 
on  the  electron,  the  potential  V  at  the  metal-oxide  interface  is  given  by 
Poisson's  relation : 


&*¥  frlTne 

(5) 

where  4  is  the  dielectric  constant  of  the  oxide.  Seaside  ring  the  term  y 
*®3Se  the  total  oxide  thickness  at  any  given  tins,  then  integrating  (5) 


we  have: 


g.iVn.,.0.  (6) 

Since  the  field  dV/dy  equals  sero  at  the  metal  surface  when  the  asm 
thickness  of  the  constant  charge  density  film  is  €,  it  follow*  thmrt 


4f  a» 
* 


V .  dim  iiy.  £ ) ;  c 


Nov  the  potential  at  the  metal  surface  ia  -(&  -  v)  for  the  condition 
that  the  oxide  film  thickness  y  approaches  aero,  where  e  is  the  energy 
required  to  remove  an  electron  from  the  metal  and  place  it  in  the  oxide 
at  the  metal-oxide  interface,  and  v  is  the  apparent  positive  potential  at 
the  metal  surface  created  by  the  growing  oxide.  Ihe  term  v,  in  general, 
is  numerically  equal  to  the  potential  difference  between  the  metal-oxide 
interface  and  the  oxide -oxygen  Interface  for  metal-oxide  systems  in  which 
lattice  diffusion  of  oxygen  ions  toward  the  metal  is  negligible.  As  will 
be  shown  later,  v  is  approximated  by  the  electron  affinity  of  oxygen 
adsorbed  on  the  oxide)  it  is,  therefore,  a  constant  independent  of  oxide 
thickness  and,  like  the  work  function,  it  is  substantially  Independent  of 
temperature.  3hen  C2  -  -(#  -  v). 

for  values  of  y  such  smaller  than  1,  the  term  y*/2  can  be  neglected. 
Hence,  from  Eq.  (7),  the  potential  affecting  escape  of  an  electron  from 
the  metal  surface  in  contact  with  an  oxide  of  increasing  uniform  negative 
space  charge  is  approximated  by: 

v  .  iy-V-  r)  (8) 

Bate  of  oxide  growth  dy/dt  is  proportional  to  rate  of  escape  of 

positive,  ions,  which  in  turn  is  controlled  by  rite  of  electron  escape. 

» 

Electron  current  across  the  contact  of  a  metal  and  a  eemi-condoetiag 
oxide  is  found  to  follow  an  equation  of  the  type^^i 


t ^"Crystal  Rectifiers",  H.  florrey  and  C.  Miltmsr,  p.  22,  80 -8b,  MeOraw- 

mu,  m  tm,  *9*8. 

. in:“‘  ' 

i  -  A*  sxp  ^  (exp  2^1  -  1)  (9) 

«mk m  A'  ic  a  eenetent,  i  is  the  work  function  of  the  gttal  modified  by 


contact  with  the  oxide,  and  v*  la  the  applied  voltage  across  the  uetal -oxide 
interface .  The  above  equation  satisfactorily  describes  the  observation  that 
log  1  for  a  rectifying  contact  is  linear  with  v*  over  a  Halted  range  of  applied 
voltage.  Since  In  our  present  considerations,  v*  is  found  to  be  In  the  order 
of  volts,  the  ten  exp  is  such  greater  than  unity  (at  1  volt, 

Sgpr  ■  39)|  m§£  Jience,  equation  (9).  simplifies  to  ,  . 

1  ■  A*  exp  *  ({■-“  V-^  (10) 

This  equation  neglects  the  reverse  electron  current  from  oxide  to  metal 
and  considers  as  important  only  current  leaving  the  netal*.  A  simplification 
of  this  kind  Is  Justified  by  the  high  positive  field  at  the  metal  surface 
when  a  metal  first  oxidizes  and  whenever  the  oxide  film  is  relatively  thin. 

Then  since  the  rate  of  film  growth  is  proportional  to  electron  flov,  fc-*. 
where  p  is  a  constant,  and  employing  Eq.  (8)  for  the  actual  potential  at  the 

II  t . 

metal  surface,  i.e.,  substituting  v  -  ■■■  7-™  1  y  for  vf,  we  haver 

®  t.- 


where  A  nov  replaces  A'p.  This  equation  Implies  that  the  effective  applied 
voltage  v  across  the  metal  oxide  surface  is  decreased  by  the  accumulating 
negative  charge  in  the  oxide  as  the  oxide  grows.  Integrating  on  the 
condition  that  when  y  *  0,  t  ■  0: 

y  -  4 exp  "e^kf"V^  ^  +  ^  C12) 


*Kq.  (10 )  applies  in  theory  only  to  electron  current  leaving  the  netal  and 
entering  the  oxide  when  this  is  the  direction  of  easy  flov  or  forward 
direction,  as  is  the  case  when  copper  oxidizes.  Oxidation  of  zinc,  to  the 
contrary,  corresponds  to  electron  flow  in  the  blocking  direction,  but  even 
here  current-voltage  data  for  rectifying  contacts  indicate  that  an  equation 
applies  of  the  same  fbrm  as  (9)  wiihsmt  alfesi'abluir  tarstpz  pf~ths~fplAed 
volte ga,  but  with  values  of  A*  and  differing  from  corresponding  values 
for  current  flov  in  the  forward  direction.  See  Bef .  52,  p.  8fc-85 
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which  Is  the 


ton  u  the  logarithmic  equation  observed  experimentally, 


y  -*0  In  (^  +  l) 


(13) 


^  *  VS  .  _  *  kT  .  eW  -  y)  k0  «(rf  -  y) 

**”*  “o  •  rrrsn  “* T  "  rifM«T'x  'Ip  t*  *  r  **p  k* 

The  above  equation  (12)  applies  only  to  a  Halted,  film  thickness  adjacent 
to  the  metal.  At  the  stage  of  oxidation  where  the  existing  interface 
electric  field  and  oxide  structure  arc  no  longer  conducive  to  filling  all 
available  sites  for  trapped  charge  in  the  oxide,  the  available  sites  become 
only  partially  occupied.  Bence,  the  excess  charge  density  in  oxide  further 
removed  from  the  metal  is  no  lopger  saturated  and  uniform  but  is  influenced 


by  the  potential  existing  at  any  point  in  the  film*  On  this  premise,  non- 
uniform  charge  density  becomes  a  property  of  films  that  are  thicker  then 
those  considered  above.  The  situation  is  similar  to  the.  diffuse  double 
layer  at  metal  surfaces  in  contact  vith  electrolytes,  as  described  by  Gouy^^, 


*53'0.  douy^  J.  de  Physique, 457  (1910). 

Chapmen^)  and  Btenr3^.  for  oxide  films  still  dd 

efeer  than  the  diffuse 

WD.  Chapman,  Phil.  Msg.,  2J,  475  (1913). 

<55>0.  Stem,  2.  naktrochem.,  JO,  506  (1924). 

barrier  laymr,  the  space  charge  disappears,  and  the  cadis  becomes  electrically 
neutral. 


When  the  charge  density  in  the  oxide  depends  on  potential  T,  the 
density  of  positive  charge  in  accord  vith  the  derivation  of  Mott  and  Quraey^^ 


can  be  expressed  by  nx  exp  $  and  the  corresponding  density  of  negative 
eV 

charge  by  n2  exp  where  ni  and  ns  respectively  represent  the  density  of 
lattice  sites  at  which  positive  and  negative  charge  nay  attach  itself.  Die 
excess  negative  charge  density  a  is  the  difference  of  these  expressions,  sad 
assuming  m.  ■  n*  ■  nQ 


n  a  nQ  {exp  -  exp  )  (lJ0 

If  the  energy  eV  is  small  compared  with  kT,  the  exponential  terms  can  be 
expanded  neglecting  higher  power  tern  then  the  first,  whereupon, 


2  n  eV 
o 

“St — 


(15) 


Substituting  this  expression  for  n  into  the  Bclsson  equation  and  assuming 
the  potential  to  be  zero  when  the  field  is  zero,  and  V  «  Ci  at  T  *  6,  where 
Y  marks  the  distance  into  the  diffuse  charge  density  layer: 


-2Y  2  n  e  Ci  _2y 

Ci  exp  and  n  »  — '■■■■"  e*P  -jr* 
o  o 


where  x  is  equal  to  and  Ci  is  the  potential  at  the  oxide- 

O  C  H  X*  C 

o 

oxygen  Interface  corresponding  to  the  oxide  thickness  L  marking  a  transi¬ 
tion  between  uniform  and  diffuse  charge  density.  Die  density  of  negative 
charge  n  obviously  falls  off  rapidly  with  distance  Y  into  the  oxide 
(Fig.  3).  If  we  now  consider  the  potential  at  the  metal -oxide  interface 
produced  by  the  diffure  space  charge,  potential  V  must  become  more  negative 
with  increased  total  thickness  of  oxide  y,  following  an  equation  identical 

j 

• 

with  the  above,  except  for  sign  of  the  exponent  and  where  Cj  is  replaced 
by  Cq  equal  to  the  potential  at  the  oxide -metal  interface  when  y  •  0.  If 
we  assume  for  simplicity  that  y  *  0  when  the  thickness  of  oxide  is  L  corres¬ 
ponding  to  a  transition  from  uniform  to  diffuse  charge  density,  then  fp am 
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Uniform  Charge! 
Density  | 

Barrier  Layer 


Diffuse  Charge  Density 
Borrler  Layer 


I 

|  Space  Charge  Loyer 

OXIDE 


H&l  2  Schematic  Diagram  of  Trapped  negative 
Cbarge  in  Utaifora  and  Diffuse  Barrier 
layers.  Positive  charge  aiy  he  localised 
as  show,  or  distributed  throughout  a 
larger  volume  of  oxide. 


*4.  (8)* 


„  *n:  / 


<=„  .  nm.Jt  .  .  t)  (M) 

*B»e  oxidation  rat*  as  before  la  proportional  to  the  nusber  of  electrona 
escaping  from  the  netal  into  the  oxide  per  unit  tills  t  Ptoftlqi  (10)  ait 
(l6),  and  the  relation »  V  *  CQ  exp  |*  ,  the  oxidation  rats  Is  tfvan 
by  the  fblicMng  exjdwseioni  °  vm  •vr^.ecwm.t 


dT  v-®?  . . . :'  ^  :v  ' 

^  ■  A  exp  — -^gy — *-  (l7) 

Furthermore,  on  the  condition  that  y  «  the  tern  exp  can  bit  expended 

*  xo 

and  approximated  by  1  +  St  .  As  will  be  shown  later,  n  for  CU4O  at  225 *C 

ls>  0  xo  • 

equals  If, 9  x  1 <r  ,  €  «=  10. 5,  and  -g-  ,  therefore,  equals  16,000  A.  This 

thickness  is  measured  not  from  the  metal  surface,  but  from  the  point  at 

which  the  oxide  grows  beyond  uniform  density  of  negative  charge;  it  Is 

appreciably  larger  than  values  of  oxide  thickness  y  with  which  we  ate 

presently  concerned,  and  cm  this  basis  expansion  of  the  exponential  term 

la  justified.  Integrating  on  the  condition  that  when  y  •  6,  t  »  0, 


x„  -C  e  C  e  - 

y  "  T  *“  [(jr  ’ll"  exp  +  j  (38) 

o  o 

which  again  is  of  logarithmic  type:  y  ■  kQ*  In  (^,  +  l)  yhafs  k^1  • 

xrt  v*r  *  Wm  -C_e  k  *  -Ci 

T  and  r*  •  ^  (>^~)  exp  -j-  exp 

0  0 

j 

• 

We  obtain,  therefore,  two  logarithmic  equations;  the  first,  Sq.  (12), 
expressing  the  thickness  of  oxide  y  as  a  function  of  tiMe  daring  build  up 
of  space  charge  having  uniform  charge  density,  and  the  second,  lq^  (18), 
corresponding  to  build  up  of  oxide  having  diffuse  space  ehmfgs.  .BtftfK 


-  ej 


e  r 
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two  equations  Indicate  that  a  change  of  oxidation  rate  should  occur  at  the 
transition  from  the  first  type  space  charge  to  that  of  the  second.  Data  In 
the  literature  support  this  conclusion!  discontinuous  changes  in  the  logarithmic 
oxidation  rate  having  been  observed  for  several  ms  tala,  including  copper^ 
iron^2'®\  and  sinc^.  Recent  data  for  copper  were  obtained  for  the 
temperature  range  150-250*0  by  MacNaim^^.  Fig.  4  is  a  plot  of  his  results 


*  *  i 


'J.  MacNaim  and  H.  H.  Uhlig,  unpublished  data. 


at  225*0 i  the 


data  at  other  temperatures  show  a  similar  linear  rela¬ 


tion  between  y  and  log  t,  and  an  increase  of  the  slope  to  a  higher  value 
at  a  critical  thickness  of  oxide.  Specimens  of  QFHC  copper  measuring  lxl 
x  3/8  inch  were  heated  in  oxygen  using  a  resistance  winding  cemented  to  one 

face.  A  thermocouple  was  embedded  in  the  metal  specimen,  both  for  measuring 

*  ■ 

temperature  and  regulating  It  automatically.  The  specimen,  with  surface 
prepared  using  a  clean  NO.  4  file,  was  introduced  into  an  all -glass  chamber 
containing  flowing  purified  Nz  and  the  temperature  brought  up  to  the 
required  value  using  low  voltage  A.C.  current.  Oxygen  was  then  admitted, 
and  the  time  of  oxidation  measured  from  this  point,  after  which  the 
specimen  was  cooled  in  Ns.  Thickness  of  the  Cu20  film  was  determined  by 
cathodic  reduction  in  0.1  N  KC1  employing  a  small  constant  D.C.  current. 

Observed  values  for  k  .  k  T  and  T"  are  summarised  in  Table  U, 
Values  of  T  were  obtained  by  extrapolating  the  first  slope  to  y  -  0,  and 
for  r  *  by  noting  the  time  at  which  change  of  slope  occurs.  In  each 
instance,  the  observed  value  of  kQ’  is  consistently  higher  than  kQ.  Values 
of  r,  inherently  subject  to  large  experimental  error,  vary  from  0.03  to 
0.3  minute  averaging  Q.12  and  for  f .«  vary  from  6  to  Xk  minutes, 
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averaging  11  minutes.  Thickness  of  oxide  L  at  which  a  change  of  slope  occurs 
and  vhich  Is  equivalent,  according  to  the  foregoing  theory,  to  thickness  of 


that  portion  of  the  harrier  layer  having  constant  charge  density,  Increases 
by  a  factor  of  ^.5  going  from  150#C  to  250*C.  Lustman  and  Mehl's^  oxida¬ 
tion  measurements  on  single  crystals  of  copper  scatter  appreciably,  but  show 
similarly  that  kQ'  is  consistently  greater  than  kQ,  with  values  generally 
of  the  same  order  of  magnitude  as  data  of  Table  II  for  overlapping  tempera¬ 
tures.  Their  values  for  T  and  Tf  vary  from  0.02  to  0.8  minute  and  from 
8  to  80  minutes,  respectively,  within  the  temperature  range  80*  to  155°C. 

The  significance  of  these  terms  will  be  discussed  shortly. 

The  Cubic  Equation 

It  Is  of  Interest  to  note  that  If  the  density  of  negative  charge  n. 

In  the  diffuse  barrier  layer  falls  off,  not  as  described  In  deriving 
Eq.  (18),  but  In  accord  with  the  Boltzmann  expression: 

n  -  nQ  exp  (19) 

then  upon  substituting  into  the  Poisson  equation,  the  solution  for  potential 
V  Is  of  the  type^'*®^: 

WRef.  45,  P*  170. 


V.lLW  +1)  (20) 

e  *0 

where  as  before  xQ  «  .  The  show  solution  assumss  that  the 


If  Instead  ws  place  the  condition  that  at  the  artel  surface  V  -  CQ  then 
y  »  0,  where  CQ  is  the  potential  at  the  nets!  surface  for  critical  Ifeiakaase 
of  unlfofk  harrier  layer  L,  then  - 


y  .  Aa  (^-  +  exp  -g-y)  (2l) 

o 

The  sign  of  the  potential  la  such  that  V  becomes  more  negative  the  greater 
the  oxide  thickness  y.  Substituting  Into  2%.  (10)  for  electron  current 
from  metal  to  oxide,  which  is  proportional  to  oxidation  rata,  wa  obtain: 

|-A«*pg  or  (22) 

■  A  exp  |n  (£-  +  exp  -gpp)  2  (2j) 

o 

Integrated  on  the  condition  that  y  -  0,  t  •  0,  we  obtain  an  equation  of  the 
cubic  form: 

(y*0«p-sr (25i) 

0 

The  cubic  oxidation  equation  baa  been  shown  to  hold  approximately 
within  certain  ranges  of  temperature  for  oxidation  of  copper 

159>T.  fihodin,  J.A.C.S.,  J2,  5102  (1950). 

^R.  Tyiecote,  J.  Inst.  Metals,  §1,  68l  (1953). 

nickel titanium^2 and  tantalum^ One  is  led  to  the  conclusion, 

Engell,  K.  Hauffe  and  B.  Ilschner,  Zelt.  llaktmoten.,  58,  478  (1954). 
^2^J.  T.  Waber,  J.  Chaa.  Ihysics,  20,  734  (1952). 

therefore,  that  charge  distribution  in  the  barrier  leper  may  follow  norm 
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than  a  single  pattern,  depending  on  experimental  conditions,  and  that  the 
particular  pattern  followed  becomes  apparent  through  the  oxidation  behavior 
of  the  metal. 

A  cubic  rate  equation  has  been  derived  by  Mott^*2^  employing  the 


'  ^'N.  Mott,  Trans.  Faraday  8oc.,  £6,  472  (1940). 

premise  that  diffusion  of  cation  vacancies  is  rate  controlling  where  the 
number  of  vacancies  is  proportional  to  rrueber  of  negative  ions  per  unit 
surface  area  and  to  a  linear  field  set  up  by  the  ions.  Engell,  Hauffe  and 
Ilschner^"^  derived  a  similar  equation  based  on  migration  of  positive  holes 
and  diffusion  of  lattice  vancancies.  The  present  derivation  is  based  on 
electron  flov  from  the  metal  as  the  controlling  step  in  the  oxidation  process, 
and,  hence,  differs  basically  from  both  derivations  of  Mott,  and  of  Bngell 
et  al.  Further  experimental  work  is  needed,  including  effect  of  grain 
orientation,  of  Curie  temperature  and  of  lattice  transformation  on  con¬ 
stants  of  the  cubic  equation  in  order  to  decide  whether  the  rate  for  any 
given  metal  Is  controlled  by  processes  of  diffusion  and  conduction  In  semi¬ 
conductors,  or  by  electron  flov  from  metal  to  oxide  at  the  stage  where  devia¬ 
tions  from  the  logarithmic  equation  become  marked.  For  thicker  films,  where 
the  parabolic  equation  takes  over,  it  is  clear  through  the  work  of  Nhgnar^ 
that  ion  diffusion  and  electric  conduction  processes  become  dominant. 
Derivation  of  the  Pideal-Jcnes  Relation 

From  Xq.  (il),  taking  the  logarithm  of  both  sides,  we  have: 

fa  -  ~~ ~  y  *  to  A  •  <*) 


Also,  from  the  relation  y  »  tn(^  +  l),  we  have: 
In  -  In  kQ  -  ln(t  +  r) 
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(8fe)  tod  (85),  and  remembering  that  kQ  ■  -  , 

kp  «  ♦  **^r+  *»  A  (*) 

Bm  oonstant  A. la  essentially  taepereturs  lnd« pendant.  Sine* 

» 

(a¥)  •  sr  ,  the  tax*  r  la  tha  time  necessary  to  oxidise  the  natal  to 

"  t  •  0 

thickness  k  at  the  Initial  oxidation  rata.  Xt  vu  found  by  Tii— ainn  an# 
K8sterv  '  to  be  strictly,  independent  of  temperature  baaed  oh  their  oxidation 
data  for  a  variety  of  aatala.  Although  it  le  coaoelvabla  that  T  nsy  actually 
varyf wltiTiaa^eratvS^lillahtly/  for  soae  metale  and  not  at  all  for  o there, 
precision  of  present-day  data,  including  values  in  Hekla  XI,  are  not  adequate 
to  Imply  a  conclusion  different  froa  that  of  Tsaeann1!.  Therefore, 

aay  he  considered  negligible  or  aero  andt 


A  fc  ko 

aT/s' 


4^1 


(27) 


vhere  AB  is  the  activation  energy  for  oxidation  in  electron  volts  froa  the 

_»AX  j 

Arrhenius  relation  kp  »•  const  exp  -gp  ,  f  is  the  work  function  of  the  metal 
modified  by  contact  with  the  oxide,  and  v  le  a  constant  positive  potential 

I 

In  the  oxide  at  the  metal  surface  ( approximated  by  the  electron  affinity  of 
the  oxygen  atoa).  If  XA  and  Xfi  are  work  functions  of  the  oxide  at  the 
Mi  oxide -oxygen  interfaces,  respectively,  it  fellow  that 


t* » -  Xx  +  ix,  -  a  (ae) 

vhere  fiQ  le  the  work  function  of  the  metal,  ^  le  the  electron  affinity 

of  oxygen  sad  ^  la  a  seel  1  additive  energy  introduced  by  adsorption  of 

a  ,  j,i  ■  •  * ' ~- 

oxygen  ea  the  oxide.  If  X,  -  X„ 

*  -  K  -  ♦  *s  )  -  *0  -  *  (89) 

vhioh  la  tin  Mdeal-Jhaee  relatiaa  aa  applied  to  enftdatiaa  of  aetale  ftamtag 


89 


w 


* 


* 


thin,  non-volatile  oxides,  and  is  the  empirical  expression  originally  proposed 
for  W,  Pt  and  C  which  form  volatile  oxides. 

If  the  work  function  of  the  oxide  at  the  metal  interface  Xi  is  less 
than  the  work  function  Xa  at  the  oxygen  interface,  then  K  will  be  larger 
than  the  value  calculated  from  Eq.  (29)  above.  This  may  be  the  situation 
for  both  Cu  and  Fe.  Values  of  i  -  AE  or  of  K  for  Cu  are  4.5  -  0.28  -  4.3 

Q 

e.v.,  and  for  Fe  are  4.4  -  0.30  -  4.1  e.v.,  where  the  activation  energy  for 
Cu  is  taken  from  MacHairn ' s  data,  end  for  iron  from  Davies,  Evans  and 
Agar's^11)  data*.  The  work  functions  are  from  Michaelson^^.  These  values 
of  K  are  larger  than  3*8  e.v.,  the  average  value  for  elements  listed  In 
Table  II.  When  more  precise  data  are  available  both  for  metal  work  functions 
and  for  activation  energies  of  oxidation,  more  definite  conclusions  will  be 
possible  than  are  now  warranted. 

Using  a  derivation  similar  to  that  given  above,  it  can  also  be  shown 


d  In  k  *  C_e 
o  o 

“TT/f — r 


(30) 


♦Since  the  logarithmic  equation  applies  to  the  oxidation  of  iron  over  the 
temperature  range  considered,  the  activation  energy  is  obtained  from  the 

slope  of  log  weight  gain  at  t  -  constant  vs  jj?.  This  follows  from  the 

proportionality  between  thickness  of  oxide  y  and  weight  gain,  and  the  fact 

that  iny«ink0+|nAn(^  +  l).  Hence,  since  T  is  a  constant  and  t 

is  held  constant,  -  H  y  -  where  k  is  Boltsmann's  constant. 

d*  dT  * 

*  fll) 

Data  in  Fig.  13  of  Davies,  Evans  and  Agar's  paper'  '  show  that  essentially 
the  same  activation  energy  is  calculated  for  various  times  of  oxidation 
(13  minutes  to  8  hours)  in  the  logarithmic  range,  with  the  exception  of  the 
13-  and  30 -minutes  runs  where  experimental  error  is  perhaps  greatest. 
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where  Cg  Is  given  by  Eq.  (l6),  and  k^*  is  the  second  sad  steeper  slope  obtained 
by  plotting  y  vs  In  t. 

DISCUSSION 

In  accord  with  Equations  (27)  and  (30),  plots  of  log  kQ  and  log  kQ*, 
using  data  of  Table  II,  vs  ^  provide  means  for  calculating  (0  -  v)  and  -CQ, 
respectively.  From  slopes  of  data  so  plotted  in  Fig.  5,  these  values  are 
0.20  e.v.  and  0.42  e.v.,  respectively.  Since  CQ  *  ~ £ji  _  (0  -  v) 

(Eq.  16),  we  knov  all  terms  of  this  equation  except  the  product  n  I,  for 
which  we  can  solve.  On  the  basis  that  &  a  10. 5  and  observed  L  at  225 °C  * 

1540  A  (Table  IB),  n  I  «  5*3  x  1010.  Also,  from  the  relation  k  -  r 

Q  »  Tv  u  C  JL 

and  employing  the  observed  value  for  kQ  at  225*0  equal  to  310  A,  we  obtain 

n  1  ■  8.1  x  1010  which  is  in  reasonable  agreement  with  the  value  calculated 

from  C  and  (rf  -  v).  We  can  next  calculate  a  value  for  n  from  the  relation 
o  o 

X  VT 

given  in  Eq.  (l8)  that  kQf  ■  -7^  ),  and  from  the  observed  value  of  kQr  ■ 

o 

o  S  k  T  & 

1630  A,  where  xq  -  (jpjfn — £2"V  •  3h*  corresponding  value  for  nQ  is  4.9  x 

o 

12 

10  ,  vhich  is  the  density  of  available  sites  for  trapped  electric  charge 
in  the  diffuse  space  charge  layer.  This  density  is  equivalent  to  about 
one  lattice  site  for  every  5  x  10^  molecules  of  Cus>0. 

We  have  next  to  check  how  density  nQ  compares  with  density  n  of 
trapped  electrons  in  the  uniform  space  charge  layer.  Since  only  the  product 
n  1  is  known  from  experiment,  an  estimate  of  1  must  be  made  through  Eq.  ,(/f) 
substituting  V  -  -ff  at  y  «  i  and  C2  ■  -(0  -  v).  It  follows  that  v  * 

i  .  Since  the  value  of  v  is  the  same  order  as  K  (Eq.  29)  or  the 
same  order  as  the  modified  electron  affinity  of  oxygen  equal  to  3.6  electron 


volts,  vs  find  substituting  6.7  x  10'iJ  for  the  average  value  of  n  I,  in 
accord  with  the  foregoing  calculations,  that  {L*  6.2  x  10  cm.  This  value 
for  1  represents  the  theoretical  maximum  thickness  of  the  uniform  charge 
density  barrier  layer  (the  diffuse  barrier  layer  may  be  considerably  thicker). 
It  is  larger  than  the  observed  thickness  of  uniform  barrier  layer  L  by  a 
factor  of  about  40,  and  should  be  approximately  comparable  with  the  thick¬ 
ness  of  the  "blocking"  layer  next  to  the  copper  surface  in  the  Cu-CugO 
rectifier.  The  thickness  of  the  blocking  layer  is  estimated  from  capa¬ 
citance  measurements  to  be  about  X0~^  cm^^,  vblefc  if  net  far  from  the 


'W'Bef.  45,  page  178. 
above  calculated  value  for  J. 


.  6.7  X  1010 


Knoving  the  value  for  1,  the  value  for  n,  therefore,  equals 


6.2  x  10' 
14 

or  1.06  x  10  .  Accordingly,  the  actual  density  of  trapped  charge  in  the 

Cu80  layer  of  thickness  L  next  to  the  metal  appears  to  be  larger  than  the 
density  (4.9  x  1012)  of  available  sites  for  trapped  charge  in  the  diffuse 
space  charge  layer  farther  away  from  the  metal  surface . 

This  relative  difference  in  density  persists,  although  absolute 
densities  increase,  if  correction  is  made  for  the  larger  absolute  surface 
compared  to  the  geometrical  surface  of  a  filed  copper  surface.  Assuming 
a  roughness  factor  of  3,  comparable  with  published  values  of  3.  8<65) 

and 


'6''K.  Powers  and  N.  Backerman,  J.  Slectrochea.  Soe.,  100,  314  (1953). 
2^(66)  tQT  abraded  metal  surfaces,  film  thickness  values  such  as  those 


'0.  Brbacher,  Zeit.  Fhysik.  Cham..  163.  215  (1933);  Chemiker  Belt., 
fc,  601  (1938). 
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in  Fig.  4. would  be  decreased  to  l/3  the  given  values,  and  both  nQ  and  n 
would  Increase  by  a  factor  of  9»  asking  the  densities  approximately  5  x  10» 
and  101-5,  respectively.  Values  of  ^  and  L  correspondingly  would  be  reduced 
by  1/3. 

Values  for  n.  i.  and  n  calculated  from  MacKairn's  data  at  several 

o 

temperatures  and  not  corrected  for  roughness  factor  are  given  in  Sable  HZ. 
Values  of  nl  are  averages  of  values  obtained  from  observed  k0  and  L,  with 
the  exception  of.  250*C  where  only  the  value  obtained  from  L  was  used.  This 
procedure  was  followed  because  the  value  for  kQ  at  250*C  falls  off  the  line 
in  Fig.  5.  If  the  value  for  kQ  were  Included,  the  corresponding  calculated 
value  for  nl  is  0.9  x  1011  compared  with  0.5  x  10^  given  in  Table  HI.  Xt 
is  evident  that  for  cuprous  oxide  in  contact  with  copper,  nQ  decreases  with 
temperature  (fewer  sites  for  trapped  charge  In  the  diffuse  barrier  layer  as 
the  temperature  Increases).  This  is  also  true  of  n,  the  actual  charge  density 
in  the  uniform  space  charge  layer.  Accordingly,  the  equilibrium: 

Bound  hole  +  free  electron  tz*  trapped  electron 
is  displaced  to  the  left  as  the  temperature  increases.  Ihe  observed  increase 
of  1  with  temperature  parallels  extension  of  the  space  charge  to  greater 
distances  into  the  oxide  the  higher  the  temperature.  This  trend  is  also 
reflected  by  Increased  observed  thickness  L  for  the  uniform  barrier  layer 
(Table  H). 

Brat  tain  calculated  the  density  of  holes  in  Cu20  as  function  of 

Brat  tain.  Bev.  Modern  Physics.  2g,  205  (l9$l)» 

temperature  from  conductivity  of  the  oxide  and  mobility  of  carriers.  Be 
found  that  a  limiting  uppsr  density,  equal  to  the  density  of  electron 
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acceptors  sinus  density  of  "acceptors  occupied  dby  electrons  fros  higher 
lying  donor  levels,  occurs  in  the  range  100*  C  or  higher  and  that  the 
limiting  density  is  iO1^  for  OFHC  copper  end  lQ^  forOhile  copper.  Obese 
values  are  reasoaably  cosparable  with  the  calculated 'density  of  charge  n 
or  density  of  defect  sites  for  bound  electrons  ng  in  Table  in.  His  cal¬ 
culations,  however,  show  that  the  density  of 'holes  does  not  decrease  with 

i  ;  '  •  '  _ _ .  . 

temperature,  but  instead  is  relatively  insensitive  to  tesqperature  in  the 

region  of  100* C.  or  higher,  and  that  below  100*  C  the  density  increases  with 
temperature. 

The  value  for  A  in  Eq.  (ll)  .can  be  calculated  from  the  relation. follow- 
ing  Eq.  (13),  T  •  -y  exp  .  At  225°C,  for  example,  T  «  3.6 

seconds,  kQ  »  310  A,  (0  -  v)  «  0.28  volt.  Therefore,  A  -  (0.86  x  10"^) 

e6»53  M  5,9  x  20"^  ca/sec.  It  will  be  remembered  that  A  ■  A’p. i$£!9.0  _ 

expresses  the  relation  between  electric  current  and  rate  of  oxide  growth, 
and  whose  value  for  density  of  CugO  «=  6,  is  1.23  x  10 cm3/ coulomb.  The 

value  for  A1  in  Be.  (lO),  accordingly,  equals  5.9.  ■  4.8  amp. /cm2. 

1.23  x  10“** 

The  average  value  for  A*  for  all  temperatures  of  Table  H  is  2.2  amp. /cm2. 

V  “Coe 

Similarly,  the  average  value  of  A*  from  the  relation  T*  ■  -y-  exp  ^y 

■Oe 

averaging  all  values  of  T*  and  kQ*  exp  -yy  is  equal  to  3*5  amp. /cm2. 
These  two  values  are  the  ease  order  of  magnitude.  These  average  values  of 
A',  if  corrected  for  a  roughness  factor  of  3,  would  be  l/3  the  given  values. 

The  ratio  -y  as  derived  from  Bq.  (12)  and  (l8)  is  equal  to  the 
k  * 

expression  exp  (-C  e/kT  -  e(ff  -  v)/kT).  tte  exponential  factor  ie 
Ko 

a  constant  at  any  given  temperature  since  -CQ  and  (ft  -  v)  are  constants 
equal  to  0.42  and  0.28  volt,  respectively.  At  225*C,  for  example,  y  - 
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exp  3.3  ■  1^3.  The  observed  ratio  is  «  167  vhich  is  of  the  same 
order  of  magnitude  .  The  average  of  calculated  ratios  for  -£r  for  all 
temperatures  is  133  corresponding  to  an  observed  average  (Table  H)  of  138 
showing  that  the  agreement  of  theory  with  experiment  is  good. 

The  described  dependence  of  oxidation  rate  on  work  function  of  the 
metal  and  on  properties  of  semi -conducting  oxides  in  contact,  immediately 
suggest  several  reasons  for  experimental  hurdles  usually  encountered  in 
reproducing  thin -film  oxidation  behavior.  Despite  refined  techniques  for 
metal  surface  preparation  and  gas  purification,  it  is  extremely  difficult 
to  avoid  the  trace  Impurities  which  greatly  affect  the  metal  work  function, 
and  also  the  number  of  sites  available  for  trapped  charge  in  oxide  films. 

The  major  Influence  of  trace  impurities  In  semi -conducting  germanium  and 
silicon  has  only  recently  been  fully  appreciated  in  connection  with  the  use 
of  these  substances  for  rectifiers  and  for  transistors.  Added  to  the  marked 
effect  of  impurities,  often  present  In  amounts  belov  common  analytical 
techniques.  Is  the  appreciable  effect  of  thermal  treatment  on  the  density 
of  lattice  imperfections  In  semi-conductors.  It  Is  little  wonder  that 
various  Investigators  haw  reported  difficulty  in  reproducing  oxidation 

(7) 

data  at  low  temperatures j  Vernon,  Akeroyd  and  Stroud' ' '  reported  such 
difficulties  in  their  study  of  sine,  Leontis  and  Hhlnes^)  in  their  study 
of  magnesium  and  Brasunas  and  TJhlig^2^  in  their  studies  of  Gr-Js  alloys. 

A  full  appreciation  of  the  various  factors  entering  thin  film  oxidation 
should  help  such  stadias  in  the  future. 
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